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ABSTRACT 

Context. The MAGIC collaboration recently reported correlated X-ray and very high-energy (VHE) gamma-ray emission from the 
gamma-ray binary LSI +61 303 during ~60% of one orbit. These observations suggest that the emission in these two bands has its 
origin in a single particle population. 

Aims. We aim at improving our understanding of the source behaviour by explaining the simultaneous X-ray and VHE data through 
a radiation model. 

Methods. We use a model based on a one zone population of relativistic leptonic particles at the position of the compact object and 
assume dominant adiabatic losses. The adiabatic cooling time scale is inferred from the X-ray fluxes. 

Results. The model can reproduce the spectra and light curves in the X-ray and VHE bands. Adiabatic losses could be the key 
ingredient to explain the X-ray and partially the VHE light curves. From the best-fit result we obtain a magnetic field of B = 0.2 G, 
a minimum luminosity budget of ~ 2 x 10 35 ergs -1 and a relatively high acceleration efficiency. In addition, our results seem to 
confirm that the GeV emission detected by Fermi does not come from the same parent particle population as the X-ray and VHE 
emission. Moreover, the Fermi spectrum poses a constraint on the hardness of the particle spectrum at lower energies. In the context 
of our scenario, more sensitive observations would allow us to constrain the inclination angle, which could determine the nature of 
the compact object. 
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1. Introduction 

LSI +61 303 is one of the few X-ray binaries (along with 
PSRB1259-63, LS5039 and probably CygnusX-1) that have 
been detected in very high-energy (VHE) gamma rays'. This 
source is a high-mass X-ray binary located at a distance of 
2.0 + 0.2 kpc (Frail & Hjellming 1991) and contains a com- 
pact object with a mass between 1 and 4 M orbiting the main 
Be star every 26.4960 + 0.0028 d (Gregory 2002) in an eccen- 
tric orbit (see Aragona et al. 2009, and Fig. 1). The main B0 Ve 
star has a mass of = (12.5 + 2.5) M , a radius of « 
10 R G and a bolometric luminosity of 10 38 erg/s (Casares et al. 
2005). Observations of persistent jet-like features in the radio 
domain at ~100 milliarcsecond (mas) scales prompted a classi- 
fication of the source as a microquasar (Massi et al. 2004; see 
also Massi & Kaufman Bernado 2009 for a radio spectral dis- 
cussion), but later observations at ~ 1-10 mas scales, cover- 
ing a whole orbital period, revealed a rotating elongated fea- 
ture that was interpreted as the interaction between a pulsar wind 
and the stellar wind (Dhawan et al. 2006). In the X-ray domain, 
LS I +61 303 shows an orbital periodicity (Paredes et al. 1997) 
with quasi-periodic outbursts in the phase range 0.4-0.8. The 
source also shows short-term flux and spectral variability on 
time scales of kiloseconds (Sidoli et al. 2006; Rea et al. 2010). 
LS I +61 303 has been detected in the VHE domain by MAGIC 
(Albert et al. 2006) and VERITAS (Acciari et al. 2008). It shows 
a periodic behaviour (Albert et al. 2009) with maxima occurring 
around phase 0.6-0.7 and non-detectable flux around periastron 



1 HESS J0632+057 is a variable galactic VHE source thought to be 
a binary, but no orbital period has been found yet (Hinton et al. 2009; 
Falcone et al. 2010). 



(0per = 0.275). However, it must be noted that since the be- 
ginning of 2008 there have been no reports of VHE detection 
around apastron even though VERITAS performed several ob- 
servation campaigns (Holder 2010). A recent detection of the 
source around periastron (Ong 2010) indicates that its VHE or- 
bital light curve might not be stable over long periods. Recently 
LS I +61 303 was detected by the high-energy (HE) gamma-ray 
instrument Fermi/LAT (Abdo et al. 2009). Its emission at GeV 
energies was found to be anti-correlated in phase with the X-ray 
and VHE emission and presented a spectrum compatible with a 
power law and an exponential cutoff at ~6 GeV. This spectrum 
is reminiscent of pulsar magnetospheric emission, but no pul- 
sations were found and no orbital variability would in principle 
be expected from this kind of emission. Since March 2009 the 
source shows a higher phase-averaged flux than before and the 
orbital modulation has significantly diminished (Dubois 2010). 
This change in GeV behaviour, as well as the change in the VHE 
light curve, indicates that a long term variability effect may be at 
play. It is not yet clear, however, whether there is any relation to 
the ~4.6 yr superorbital modulation of the radio peak (Paredes 
1987; Gregory 1999). 

Two scenarios have been proposed to account for the ob- 
servational features of LS I +61 303. The microquasar scenario 
requires accretion onto the compact object and the formation 
of a jet in which particles are accelerated and emit up to 
VHE gamma-rays through upscattering of stellar photons (e.g. 
Bosch-Ramon & Paredes 2004; Bosch-Ramon et al. 2006a). In 
the pulsar wind zone scenario, the compact object is a young 
pulsar and particles are accelerated at the interaction zone be- 
tween the pulsar wind and the stellar wind (Maraschi & Treves 
1981; Dubus 2006). The stellar wind is thought to be formed by 
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a fast, low-density polar wind and a slow, high-density equato- 
rial decretion disc, which may be strongly truncated and have its 
properties disturbed by the compact object (Romero et al. 2007). 
On the other hand, the characteristics of the interaction of pulsar 
winds or microquasar jets with a massive star wind are not well 
constrained (but see the studies by, e.g, Bogovalov et al. 2008; 
Perucho & Bosch-Ramon 2008). Doppler boosting effects could 
be present both for a microquasar jet and a shocked pulsar wind 
outflow (e.g. Khangulyan et al. 2008a; Dubus et al. 2010), but its 
effects in both scenarios are difficult to assess because of our lack 
of knowledge on the geometry and orientation of the emitter. 
Finally, the uncertainty in the orbital geometry, in particular the 
inclination, also affects the ability of models to recover the phys- 
ical parameters of the emitter (Sierpowska-Bartosik & Torres 
2009). All this shows the difficulties of the study of gamma-ray 
binaries. 

Because of the many factors involved, simultaneous X- 
ray/VHE observations are a useful mean to start to disentangle 
the physics of the high-energy emission in gamma-ray binaries. 
However, the combined effect of short-term variability in the X- 
ray domain and night-to-night variability in the VHE domain has 
precluded a clear detection of X-ray/VHE emission correlation 
from archival observations. In 2007, a campaign of simultane- 
ous observations with the MAGIC Cherenkov telescope and the 
XMM-Newton and Swift X-ray satellites revealed a correlation 
between the X-ray and VHE bands (Anderhub et al. 2009). The 
suggestion that the emission in both energy bands comes from 
the same population of accelerated particles (Anderhub et al. 
2009, see discussion) makes these observations an ideal data set 
to test the properties of the X-ray/VHE emitter in LS I +61 303. 
It is worth noting that the need of only one leptonic population 
to explain the X-ray and the VHE emission would be against 
lepto-hadronic models (e.g. Chernyakova et al. 2006). 

In this work, we present the results of the application of a 
leptonic one zone model to explain the measured light curves 
and spectra from these observations. Our goal is to obtain the 
physical parameters of the emitter robustly by making as few 
assumptions as possible. 

2. Simultaneous X-ray/VHE observations 

LS I +61 303 was observed during ~60% of an orbital period si- 
multaneously in the VHE and X-ray bands in September 2007 
(Anderhub et al. 2009). During a first part of the campaign 
(0.43 < 4> < 0.7), the MAGIC Cherenkov telescope observed 
the source for three hours every night and simultaneous XMM- 
Newton observations were performed. During a second part 
(0.71 < (f> < 1.13) the X-ray observations were performed with 
Swift/XRT. The observation times were of about 15 ks for the 
XMM-Newton observations and of about 3 ks for the Swift/XRT 
observations. The longer observations and higher effective area 
of XMM-Newton resulted in much higher statistics than the mea- 
surements from SwiftPiRT. To take into account the variability 
that the source shows on short time scales in the X-ray band 
when comparing these fluxes to the VHE measurements, the rms 
count rate variability of the X-ray light curve was considered in 
addition to the statistical uncertainty in the flux measurement 
to obtain realistic flux uncertainties (see Anderhub et al. 2009, 
for details). Using this method, the X-ray flux uncertainties take 
values between 6% and 25% of the flux. The simultaneous cam- 
paign revealed very similar light curves in the X-ray and VHE 
bands: a first outburst at </> = 0.62 with a rise time below 20 h and 
a decay of about two days; and a second broader outburst span- 
ning a few days between phases 0.8 and 1 . 1 with lower peak flux 
than the first. The correlation coefficient between the two bands 




Fig. 1. Orbital geometry of LS I +61 303 looking down on the 
orbital plane showing the orbit of the compact object around the 
Be star in units of the orbital semi-major axis. The positions la- 
beled with the orbital phase correspond to periastron, apastron, 
and inferior and superior conjunctions. The X-ray/VHE multi- 
wavelength campaign covered the phase range shown in green 
(0.43 < (f> < 1.13), while the phase range of the VHE emission 
outburst is shown in red (0.6 < <p < 0.7). The Be star (blue cir- 
cle) and the equatorial wind (grey circle) are shown assuming an 
equatorial disk radius of = 7R*. 



for the first outburst was found to be r — 0.97, whereas taking 
all simultaneous observations lowered the value to r — 0.81. 

Although all X-ray observations resulted in a clear detec- 
tion, some of the VHE flux measurements have a significance 
below la. For these cases we will consider the 95% CL (con- 
fidence level) upper limits calculated by Jogler (2009). The de- 
tailed fluxes and uncertainties of the X-ray and VHE light curves 
may be found in Tables 1 and 2 of Anderhub et al. (2009) and the 
VHE upper limits in Table A.l of Jogler (2009). 

The sensitivity of VHE observations precludes obtaining 
nightly spectra, and only a spectrum combining the obser- 
vations with phases between 0.6 and 0.7 was published by 
Anderhub et al. (2009). In order to obtain a simultaneous SED, 
we extracted an XMM-Newton spectrum of the source for the 
three X-ray observations that correspond to the MAGIC spec- 
trum. We filtered the data using SAS vlO.O and extracted three 
individual spectra from the pn instrument. We used the SAS tool 
efluxer to convert the spectra from counts vs. channel into 
physical units (flux density vs. energy), averaged the three un- 
binned spectra and then grouped the bins to a signal-to-noise 
ratio of 20. In order to compare the measured spectrum to the 
computed intrinsic spectrum, we deabsorbed the former using 
the wabs FORTRAN subroutine provided with Xspec vl2.0 
( Arnaud 1996). We used a column density of Nh = 5 X 10 21 cm -2 
obtained from the average of the individual fits of an absorbed 
power law to the three observations. This column density is con- 
sistent with that of the ISM alone (Paredes et al. 2007). 
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3. Model description 

The discovery of a correlation between the X-ray and VHE 
bands points towards a common mechanism of emission mod- 
ulation at both bands. In a leptonic scenario, the fast and simul- 
taneous changes in flux in both bands indicate that the mod- 
ulation mechanism has to directly affect the emission level of 
the Inverse Compton (IC) and synchrotron processes. Assuming 
constant injection, a way to obtain correlated X-ray and VHE 
emission is through a modulation of the number of emitting par- 
ticles by dominant adiabatic losses, which would be ultimately 
related to the (magneto)hydrodynamical processes in the accel- 
erator and emitter regions. These processes may be related, for 
instance, to the interaction of the pulsar wind or the black hole 
jet with the stellar wind of the massive companion. The hard 
X-ray spectrum with photon index r* 1.5 also points towards 
dominant adiabatic losses, which imply an injection electron in- 
dex of a e ~ 2. In the regime of dominant adiabatic losses the 
emitted X-ray flux is proportional to the number of emitting par- 
ticles, so the orbital dependency of adiabatic losses can be in- 
ferred from the X-ray light curve. Dominant adiabatic losses are 
not uncommon when modelling gamma-ray binaries. They have 
also been invoked to explain the variations of the X-ray and VHE 
fluxes in the gamma-ray binaries PSR B 1259-63 and LS 5039 
by Khangulyan et al. (2007) (see also Kerschhaggl 2011) and 
Takahashi et al. (2009), respectively. We note that, as explained 
in Sec. 3.4, the IC emission will be additionally modulated by 
the changes in seed photon density along the orbit as well as the 
geometrical effects of the location of the emitter along the orbit. 

For the sake of simplicity, and given how little we know 
about the exact properties of the particle accelerator and the non- 
thermal emitter in LS I +61 303, we will here assume a constant 
particle injection spectrum and emitter magnetic field along the 
orbit. We will also adopt one leptonic population cooling down 
under variable adiabatic losses at the location of the emitter, 
which will be approximated as homogeneous and point-like. 

The orbital parameters were adopted from Aragona et al. 
(2009) with an inclination angle of i = 45°. We will also discuss 
other possibilities regarding the inclination angle, in particular 
the extremes of the allowed range 15° S5 i S5 60° considered 
by Casares et al. (2005). A sketch of the orbital configuration of 
LS I +61 303 can be seen in Fig. 1 with the phase ranges of the 
simultaneous X-ray/VHE campaign indicated. The equatorial 
wind disk in binaries of the characteristics of LS I +61 303 is ex- 
pected to be truncated at a radius of rj = 5R+ (Grundstrom et al. 
2007). However, the variability in the equivalent width of the Ho- 
line indicates that its radius can vary between 4.5/?+ and 7/?+ 
(Zdziarski et al. 2010; Sierpowska-Bartosik & Torres 2009). A 
variable red shoulder in the Ha line (McSwain et al. 2010) could 
point towards the creation of a tidal stream in the equatorial disk 
because of the proximity of the compact at periastron that even- 
tually falls back onto the disk at later phases. Even at its largest 
possible radius of 7R*, as assumed in Fig. 1, the equatorial disk 
does not affect the phase ranges of the observations we are con- 
sidering, and thus we will not include it in our model. 

3. 1 . Adiabatic cooling 

One of the remarkable features of the VHE light curve of 
LS I +61 303 is the lack of detectable emission during periastron 
(Albert et al. 2009). IC emission is very effective at this phase 
owing to the high stellar photon density, so the lack of emis- 
sion means that either it is strongly absorbed or that the num- 
ber of accelerated electrons falls drastically. However, the an- 
gular dependency of the pair production process, oc ( 1 - cos ijj) 
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Fig. 2. Cooling times along the orbit for adiabatic (solid), syn- 
chrotron (dotted) and IC (dashed). The synchrotron and IC cool- 
ing times are shown at E e = 10 12 eV. This energy corresponds to 
electrons that emit in the X-ray band through synchrotron and in 
the VHE band through IC. 

for the interaction rate and oc (1 - cosi/O -1 for the gamma-ray 
energy threshold (where tf> is the interaction angle of the in- 
coming photons), implies that absorption is very low for phases 
immediately after periastron even under the dense seed photon 
field of the star (see Sec. 3.5 for details on the calculation of 
pair production absorption). Inverse Compton emission also has 
a dependency on the angle between the seed photon and the 
accelerated electron and would be lower during these phases, 
but this effect is diminished at higher energies when scatter- 
ing takes place in the deep Klein-Nishina regime (see Fig. 6 in 
Khangulyan et al. 2008b). A reduction in the number of accel- 
erated particles is required to explain the lack of VHE emission 
around periastron. Although the X-ray flux during periastron is 
around half of the peak flux during the periodic outbursts, the 
non-detection of VHE emission during periastron places an up- 
per limit at only a 10% of the outburst peak VHE flux. Taking all 
this into account, we will assume that in the context of one elec- 
tron population the X-ray/VHE correlation is only valid for the 
X-ray excess flux above a certain quiescent or pedestal emission. 
Therefore, the number of emitting particles responsible for this 
excess X-ray flux is lower and their VHE emission remains at 
non-detectable levels around periastron. This reasoning is sup- 
ported by the X-ray/VHE correlation found by Anderhub et al. 
(2009), in which the constant term of the linear correlation is 
(12.2^q) x 10 _12 ergcm _2 s _1 . Considering that some of the 
X-ray observations are below this constant term, we took the 
pedestal flux value to be F pe d = 1 1 .5 x 10~ 12 erg cirT 2 s _1 . 

As mentioned above, a scenario with dominant adiabatic 
losses is a very good candidate mechanism to explain the cor- 
relation characteristics of the source. When considering a sys- 
tem with constant magnetic field and constant injected electron 
spectrum, the only modulation in synchrotron radiation is due to 
the modulation of the dominant cooling process that results in a 
modulation of the (evolved) stationary electron distribution. For 
dominant adiabatic losses and a constant injection rate, the X- 
ray flux will be proportional to the adiabatic cooling time f al j. A 
consistent calculation of f a( j would require knowledge of the ex- 
act nature of the source and the (magneto)hydrodynamical pro- 
cesses that take place there. One can however infer ? a[ j along 
the orbit by taking it proportional to the excess flux above the 
pedestal Fx - F pe d, which is independent of the electron energy. 
For the phases covered by the X-ray observations we derived a 
smooth curve following the behaviour of the X-ray fluxes. For 
the phases outside the campaign coverage we assumed a flat 
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light curve consistent with the first four observations at phases 
0.43 to 0.59 and the last observation at phase 0.13. Choosing 
the ratio f a d/(^x - ^ped) as small as possible while keeping the 
adiabatic losses dominant allows us to put a lower limit on the in- 
jected luminosity budget. This choice of the ratio f a( j/(fx _ -Fped) 
yields adiabatic loss time scales of between 35 and 500 seconds. 
Figure 2 shows the orbital evolution of the cooling time scales in- 
ferred for adiabatic and calculated for IC and synchrotron losses. 
The adiabatic cooling time is independent of the electron energy, 
but for synchrotron and IC the cooling time scales are given at 
£ e = 10 12 eV. This energy corresponds to electrons that emit 
in the X-ray band through synchrotron and in the VHE band 
through IC. 



> 



1*5 




0.4 0.6 0.8 

Orbital Phase 



1.2 



3.2. Particle energy distribution 

The radiation emitted by a population of non-thermal parti- 
cles arises from the contribution of electrons of different ages. 
As electrons loose energy through radiative or adiabatic cool- 
ing processes, an evolved particle energy distribution, differ- 
ent from the injected particle spectrum (Q(y), taken here con- 
stant in time), arises. The evolved particle spectrum n(y e ,t) 
can be obtained from Q{y) through the following equation 
(Ginzburg & Syrovatskii 1964): 



dn(y e ,t) dyn(y e ,t) 

+ : = 6(7e). 



dt 



(1) 



where y is the sum of the energy loss rates for IC, synchrotron 
and adiabatic processes and Q(y e ) is the acceleration rate. Since 
little is known about the acceleration process responsible for 
particle acceleration in LSI +61 303, we assume a constant 
particle injection spectrum along the orbit. We adopt a phe- 
nomenological electron injection spectrum Q{y e ) consisting of 
a power law function with an exponential high energy cutoff: 
Q(je) = Go Te" c exp(-y e '«c 2 /£e,max)- The maximum electron 
energy Z? e ,max is obtained from the balance of acceleration and 
energy loss rates, thus establishing a cutoff in the accelerated 
particle energy spectrum (see Sec. 3.3); a e is typically ~ 2 in 
the context of non-relativistic first order Fermi acceleration (see, 
e.g., Drury 1983); and Qq is the normalization of the function, 
which we keep constant along the orbit. 

For LSI +61 303 the adiabatic and radiative cooling time 
scales of high-energy electrons are much shorter than the orbital 
period, and thus we can calculate the steady-state electron en- 
ergy distribution along the orbit omitting the first term in Eq. 1 
and solving for n(y e ): 



ly(0)l Jy c 



(2) 



The resulting n(y e , (p) is the evolved particle spectrum for which 
the emitted radiation properties are calculated along the orbit. 

3.3. Maximum electron energy 

The VHE gamma-ray energy spectrum of LS I +61 303 extends 
up to E y ~ 10 TeV (see Acciari et al. 2009) without a break in 
the spectrum, indicating that electrons must be accelerated up 
to energies beyond 10 TeV. The acceleration of electrons has to 
be faster than the cooling times of the different cooling mecha- 
nisms so that the electrons can reach these very high energies. 
The maximum electron energy will be that for which f acc = t/ coo i. 
The acceleration time of electrons can be characterized as 



Fig. 3. Maximum electron energy £ e n 
77 = 10 (see text for details). 



along the orbit taking 



where = EjeB is the Larmor radius, Be is the magnetic field 
in Gauss and £jev is the electron energy in TeV. The dimen- 
sionless constant 77 > 1 indicates the efficiency of the acceler- 
ation process, where r\ — > 1 corresponds to the maximum pos- 
sible acceleration rate allowed by classical electrodynamics and 
77 = 2n(v s /cy 2 for parallel non-relativistic diffusive shocks in 
the Bohm regime (Drury 1983). 

If electron energy losses are dominated by adiabatic cooling, 
the maximum energy that electrons will reach is given by f acc = 
f a d, which results in 



9B G ^ _1 TeV . 



(4) 



However, at some points along the orbit in which adiabatic losses 
are smaller, synchrotron losses may dominate at high energies, 
since f syn oc E~ l whereas the adiabatic cooling time is constant 
with electron energy. In these cases, with a synchrotron cooling 
time described by 



f syn ^400£ Te ' v B G 2 s, 



the balance f acc = f syn results in a maximum energy of 



^e,m a x ~ 60 Bq 



-'V^TeV. 



(5) 



(6) 



j]r L /c « 0.1l77£ , T evB G S ' 



(3) 



Inverse Compton energy losses at the relevant energies are 
smaller than to the synchrotron and adiabatic losses. At each 
point along the orbit, all relevant time scales are evaluated and 
the corresponding £ e niax is obtained, which in this work charac- 
terizes the position of the cutoff in the injected particle spectrum 
through an exponential cutoff. The adiabatic cooling times in- 
ferred from the X-ray light curve range from a few tens to a 
few hundred seconds. When considering the requirement that 
acceleration takes place up to electron energies of 10 TeV, we 
can see that the accelerator required must be quite efficient with 
77 « 7-130 depending on the orbital phase. Since we are assum- 
ing that the injection is constant, we also expect the acceleration 
rate to be constant and have chosen a value of 77 = 10 to perform 
the energy cutoff calculations along the orbit, which would cor- 
respond to acceleration in a parallel diffusive shock with velocity 
of ~0.8c (therefore favoring a mildly relativistic shock). The or- 
bital dependency of the maximum electron energy for 77 = 10 
is shown in Fig. 3. During periods of low X-ray flux E^ m - dx has 
a behaviour proportional to f a d (see Fig. 2). During the outburst 
peaks, on the other hand, high energy losses are dominated by 
synchrotron and the cutoff is constant (following Eq. 6 for con- 
stant magnetic field). 
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3.4. Synchrotron and IC emission 

For each position along the orbit, we calculate the emitted syn- 
chrotron spectrum from the evolved particle distribution and the 
magnetic field. Synchrotron emission in the X-ray range will be 
generally characterized by a photon index of Tx = (a e + l)/2 ~ 

1.5. The orbital dependence of the electron maximum energy 
will affect the shape of the spectrum, typically hardening (soften- 
ing) it for high (low) fluxes, which correspond to lower (higher) 
adiabatic losses and therefore a higher (lower) maximum elec- 
tron energy. However, since the cutoff in the particle distribution 
is located at energies significantly higher than those responsible 
for X-ray emission, this effect is small in the X-ray energy band 
(Ar x ~ 0.2). 

We calculate the IC component of the spectrum using the 
anisotropic cross section and the stellar radiation field (black- 
body radiation at kT = 2 eV and = 10 3S erg s -1 ) as the source 
of seed photons for the interaction. Following Eqs. 20 and 21 of 
Aharonian & Atoyan (1981), the photon spectrum emitted by a 
population of isotropically distributed electrons can be described 
with a precision better than 3% (for e » eo and y » 1) by 



d 2 N(6, e) 
<9edQ 



3 f «£„ 

- — c °t eo,6>,y)de , 

16n J^Mfi) f or 



(7) 



where e and eo are the energies of the scattered and incident pho- 
tons, respectively, in units of m e c 2 , is the scattering angle, cr T 
is the Thomson cross section, « £o is the number density of the 
stellar photon field, 



and 



2(1 -cos0)y 2 [l -e/yY 



f(e,eo,6,y) = l + 



2z 2z 2 
+ 



2(1 -z) b„(l-z) b 2 g (l-z) 2 



(8) 



(9) 



where bg — 2(1 - cos6>)eoy ar, d z — e/y. For a blackbody seed 
photon distribution, the shape of the resulting emitted spectrum 
only depends on the parameter bg and the relativistic electron 
energy distribution. 

3.5. Pair production absorption 

The intense radiation field can also absorb very high en- 
ergy y-ray emission through pair production yy — > e + e~ 
(Gould & Schreder 1967). An exploration of the effects that this 
absorption can have on the observed spectrum of gamma-ray bi- 
naries can be found in Bosch-Ramon & Khangulyan (2009) and 
references therein. The differential opacity for an emitted y-ray 
of energy e is given by 



ATyy = (1 - COS ^(TyyU^ Al dC( 



(10) 



where / is the distance along the line of sight, eo is the energy of 
the stellar photon, ij/ is the interaction angle and the absorption 
cross section can be represented in the form: 



(Ty 



-P 2 ) 



(id 



where /3 = (1- l/s) 1 ^ 2 , and s = eeo(l -cos 00/2. Pair production 
can only occur for s > 1, when the centre of mass energy of 
the incoming photons is sufficiently high to create an electron 
positron pair. The cross section maximum takes place for s 
3.5, with a value of cr yy 0.2<tt, and then decreases for s » 1. 



The optical depth owing to pair production T yy is calculated 
at each point along the orbit by integrating the differential opac- 
ity (Eq. 10) along the line of sight and over the blackbody seed 
photon distribution. The attenuation factor k = exp(-r yy ) is ap- 
plied to the intrinsic IC spectrum (dashed line in Fig. 4) to ob- 
tain the spectrum emitted out of the binary system (solid line in 
Fig. 4). 



4. Results 

As seen in Fig. 6, we were able to reproduce the X-ray and VHE 
light curves obtained during the 2007 multiwavelength cam- 
paign. Although the good agreement with the X-ray light curve 
is expected because we derived the dominant adiabatic losses 
from the observed X-ray fluxes, the VHE light curve is a pre- 
diction of the model taking into account the binary geometry, 
the stellar photon field, and the non-thermal particle distribution. 
Furthermore, the spectral indices in the X-ray and VHE energy 
bands for the outburst at phases 0.6 < <f> < 0-7 are well repro- 
duced, as can be seen in Fig. 4. 

4.1. Electron population 

We found that the best fit to the X-ray and VHE photon indices 
can be obtained by taking a particle energy distribution with 
slope a e = 2.1. As mentioned above, the photon index in the 
X-ray band directly depends on the particle index of the parent 
electron population. In the X-ray band, the observations show 
a photon index between 1.58 + 0.02 and 1.66 ± 0.02 (plus a 
1 .85 ± 0.02 outlier) for the XMM-Newton observations, which is 
anti-correlated with the X-ray flux. The time-dependent energy 
cutoff given by the variability of cooling time scales provides 
a small time dependent variation of the photon index through a 
modification of the shape of the spectrum. Regarding the VHE 
observations, the measured photon index is 2.7 ± 0.3 sta t ± 0.2 sys 
for the three observations during the first outburst. Because ab- 
sorption is low at the orbital phases of the maximum, the photon 
index is determined by a e and the IC interaction angle. The par- 
ent particle population with a e = 2.1 results in computed X-ray 
photon indices in the range 1.55-1.67 and reproduces the ob- 
served anti-correlation between Fx and Fx- In the VHE band 
the adopted a e results in photon indices in the range 2.6-2.8 for 
high-flux phases, whereas the spectrum is much softer for low 
flux phases because of the lower energy cutoff. Figure 4 presents 
the SED averaged over the phase ranges of the three simultane- 
ous observations during the first outburst (phases 0.62, 0.66 and 
0.70), and shows that both the flux levels and the photon indices 
at X-ray and VHE are well reproduced. 

We found that the phase-averaged SED calculated using a 
power law with a e = 2.1 and a high-energy cutoff produces a 
flux too high in energies between 10 and lOOGeV when com- 
pared to the phase-averaged spectrum measured by FermijLAX 
(Abdo et al. 2009). The GeV spectrum and fluxes along the or- 
bit appear incompatible with one leptonic population. A way 
to avoid this excess flux is to consider the electron injection 
spectrum described by a broken power law with a harder par- 
ticle index below £ e = 4 x 10 11 eV. We found that the softest 
particle index for which the computed GeV spectrum is com- 
patible with the FermifLAT measurement is a e % 1.8. It must 
be noted that this change in the particle index at low energies 
does not affect the spectrum of particles that emit in the X-ray 
and VHE bands through synchrotron and IC, respectively. The 
FermifLAT spectrum along with the computed phase averaged 
SED is shown in Fig. 5. The phase-averaged FermifLAT can- 
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Fig. 4. Computed SED averaged over the three observation pe- 
riods during the first outburst (phases 0.62, 0.66 and 0.70), 
with synchrotron (dot-dashed), IC (solid) and non-absorbed IC 
(dashed) components. The calculations were performed with 
i = 45° and B = 0.22 G. The crosses show the XMM-Newton 
EPIC-pn spectrum averaged over three observations and deab- 
sorbed taking Nh = 5 x 10 21 cirT 2 . To compare the measured 
flux to the excess computed X-ray flux, we subtracted a pedestal 
power law spectrum corresponding to a 0.3 - 10 keV flux of 
11.5 x 10 _12 ergcm _2 s _1 . The MAGIC simultaneous spectrum 
is shown as a red bow-tie. 




Fig. 5. Computed SED averaged over the whole orbit with syn- 
chrotron (dot-dashed), IC (solid) and non-absorbed IC (dashed) 
components, calculated with the same parameters as in Fig. 4. 
The crosses and upper limits indicate the phase-averaged 
FermifLAT spectrum. The MAGIC VHE spectrum during the 
outburst is shown as a bow-tie. It must be noted that the MAGIC 
spectrum corresponds to different phases than the computed 
SED and is shown as a reference to Fig. 4. 



not be easily compared with the peak spectrum shown in Fig. 4 
because the X-ray/VHE and GeV measurements are not simul- 
taneous. Furthermore, the orbital variability of the GeV spec- 
trum at higher energies (i.e., the shape of the cutoff) has not been 
clearly established by Abdo et al. (2009), and a phase-resolved 
spectrum cannot be derived from a phase averaged-spectrum. 



4.2. Magnetic field 

The X-ray/VHE flux ratio is a very sensitive indicator of the 
magnetic field in the emitter. We found the magnetic field that 



best describes the observed light curves by performing a ex- 
ploration of the parameter space of B and the normalization 
of the injected particle spectrum Qq. For each (B, Qo) pair we 
calculated the X-ray and VHE light curves and estimated the 
goodness-of-fit to the observed flux points using a x 1 test- The 
observed light curves are best described using an ambient mag- 
netic field of B — 0.22 G. It is not trivial to set a formal un- 
certainty range on this value because of the unknown number of 
constraints that should be considered (e.g., the shape used for 
the adiabatic cooling curve). However, we see that for variations 
of the order of 0.05 G around the best-fit value of B — 0.22 G 
the observed X-ray and VHE light curves are no longer well de- 
scribed by the computed ones. 

4.3. Injection energy budget 

The adopted luminosity in injected accelerated particles is 
around 2 x 10 35 ergs _1 along the whole orbit. The computed 
emission bolometric luminosity, on the other hand, varies be- 
tween ~10 34 and ~1.5 x 10 35 ergs~'. Note that the adopted in- 
jected luminosity is only a lower limit because we arbitrarily 
chose the lowest ratio between the X-ray flux and the adiabatic 
cooling time scale. For higher adiabatic losses (while still pro- 
portional to X-ray emission), the emitted spectrum would be 
equivalent but the required injected luminosity would be higher. 

The dominant energy band in the broad-band SED of 
LSI +61 303 is the MeV-GeV region. Using the gamma-ray 
spectrum measured by Fermz/LAT (Abdo et al. 2009) we can 
obtain a lower limit on the power available for particle accel- 
eration. The observed energy flux above 100 MeV is Gioo — 
4.1 x 10~ 10 erg cirT 2 s _1 , resulting in an observed luminosity of 
2 x 10 35 erg s _1 . Given this high luminosity, it is natural to think 
than the acceleration mechanism at play is able to provide a sim- 
ilar power to the X-ray and VHE emitting electrons. 



5. Discussion 

5.1. Adiabatic time scales 

The orbital variability and range of adiabatic losses may be 
used to pose constraints on the hydrodynamical properties of the 
various scenarios proposed for LSI +61 303. Hydrodynamical 
calculations have been performed in the binary pulsar scenario 
(Bogovalov et al. 2008; Khangulyan et al. 2008a) and for the in- 
teraction of a microquasar jet with the wind of a young star 
(Perucho & Bosch-Ramon 2008; Perucho et al. 2010). However, 
for LS I +61 303 it is difficult to calculate the adiabatic losses in 
each of the mentioned scenarios because of the uncertainties in 
the properties of the Be stellar wind and the relativistic outflows. 

In our phenomenological framework, we find that adiabatic 
loss time scales range between 35 and 500 s, as seen in Fig. 2. 
Whereas it is difficult to relate these time scales to those arising 
from complex geometries (e.g. shocked pulsar wind or a micro- 
quasar jet), we can use the simple model of an expanding sphere 
to obtain insight on relevant emitter properties. The adiabatic 
loss time scale r a( j is related to the typical scale of the emitter R &m 
and its expansion velocity v exp through f a( j = ^? e m/v ex p- While the 
size of the emitting region is unknown, the fact that a one-zone 
emitter located at the compact object position can successfully 
explain the X-ray/VHE emission indicates that the region must 
be small compared with the orbital separation. Considering an 
emitting region with a radius 10% of the orbital separation we 
obtain a range of radii between ~ 3 x 10 11 cm (at periastron) and 
~ 9 x 10 11 cm (at apastron). For these sizes, the maximum ex- 
pansion velocities would be required for the shortest adiabatic 
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Fig. 6. Top: Computed VHE light curve 
(red line) and observed VHE light curve 
by MAGIC in units of 10~ 12 ph/cm 2 /s. 
Observations with a significance above 2cr 
are shown in filled squares, while 95% CL 
upper limits are shown otherwise. Bottom: 
Computed X-ray light curve (red line) and X- 
ray light curve observed during the multiwave- 
length campaign with XMM-Newton (filled 
circles) and Swift/XRT (open circles) in units 
of 10~ 12 erg/cm 2 /s. The pedestal flux is indi- 
cated by a horizontal dashed line. All error 
bars correspond to lcr uncertainties. The same 
parameters as in Fig. 4 are used for both pan- 
els. 



cooling time scale, i.e. 35 s. At apastron the expansion would 
need to be highly relativistic, with a velocity of 0.85c (similar to 
the relativistic sound speed cf V3), corresponding to a Lorentz 
factor of T = 1 .89, whereas for periastron it would be mildly 
relativistic at 0.29c. For the longer adiabatic loss time scales of 
hundreds of seconds, the required expansion velocities would be 
much lower, of the order of 0.01c-0.06c. A change in the size of 
the emitter could be caused by the change in external pressure 
from the stellar wind, resulting in a larger emitter at apastron 
(Takahashi et al. 2009). This would explain the location of the 
X-ray/VHE outbursts around apastron, but it is not clear which 
magnetohydrodynamical process is responsible for the variation 
of the adiabatic loss time scale of more than an order of magni- 
tude in 20 h. Wind clumping can be considered as a source of 
pressure variability, as suggested by Zdziarski et al. (2010), but 
it would not explain the constant location of the X-ray/VHE out- 
burst (0.6 < <p < 0.7) orbit after orbit. On the other hand, a tran- 
sition of the emitter from a region dominated by polar wind to a 
region dominated by equatorial wind, or vice versa, could give 
rise to such a fast change in wind pressure. In addition, there 
could be other yet unknown properties of the stellar wind, as 
shown by the recent detection of a variable red shoulder in the 
Ha line by McSwain et al. (2010), which could be related to a 
tidal stream in the equatorial wind falling back onto the disk. 

5.2. Energetics 

The energetics of the system, regardless of the acceleration 
mechanism, are able to provide the required injected power used 
in our model given the high flux in the MeV-GeV band (see 
Sec. 4.3). 

The pulsar-like spectrum of LSI +61 303 in the GeV band 
(power law with exponential cutoff at a few GeV) prompted 
Abdo et al. (2009) to consider a magnetospheric origin for the 
HE gamma-ray emission. In this scenario we can gain some in- 
sight on the spin-down luminosity of the putative pulsar from the 
HE gamma-ray luminosity. The observed energy flux Gioo from 
magnetospheric emission can be related to the emitted gamma- 
ray luminosity L y as 

L y = 4nfa(a,&)d 2 G m , (12) 

where fa is a beaming correction factor that depends on the ge- 
ometry of the emitter, the magnetic angle a and the Earth view- 
ing angle The emitted gamma-ray luminosity is related to 



the pulsar spin-down luminosity through the gamma-ray effi- 
ciency r\ y = L 7 /E. A trend was found from EGRET observa- 
tions that the efficiency scaled as E~^ 2 (Thompson et al. 1999), 
and this was later confirmed (at least for E > 10 34 ergs -1 ) 
by FermifLAT arriving to T] y 0.034(£/10 36 ergs" 1 )- 1 ' 2 
(Abdo et al. 2010). Using this relation, E can be inferred from 
Ly as E = L 2 /(1.156 x lO^ergs" 1 ). From the LSI +61 303 
phase-averaged gamma-ray spectrum measured by Abdo et al. 
(2009) we obtain Gioo - 4.1 x 10 _10 ergcnr 2 s _1 , from which 
L y 2 x 10 35 /n erg s _1 . Following these considerations, the 
spin-down luminosity inferred from the FermifLAT observa- 
tions is of £ =i 3.3 x 10 37 /^ erg s _1 . Traditionally (e.g. polar 
cap models) it was assumed that the gamma-ray beam covers 
a solid angle of 1 sr uniformly, so fa = l/(4n) 0.08, which 
would result in a spin-down luminosity of E 2 x 10 35 erg s" 1 . 
However, more modern pulsar beaming models (e.g. slot gap 
and outer gap models) have fn~l for many viewing angles 
(Watters et al. 2009), which would result in a more energetic 
pulsar with E 3 x 10 37 erg s . The fact that the cutoff of the 
spectrum is located at high energies (Ecutoff = 6.3 +1.1 GeV) 
points to a high-altitude location for gamma-ray emission, since 
emission beyond a few GeV is very difficult close to the neu- 
tron star because of strong attenuation from y - B — > e + e~ 
absorption (Baring 2004). The general trend in the FermifLAT 
pulsar sample (Abdoetal. 2010) also points towards a high- 
altitude location, so we favour slot gap or outer gap models 
for which the spin-down luminosity of the putative pulsar in 
LSI +61 303 would be of £ * 3 x 10 37 ergs _1 . Even though 
the FermifLAT observations do not completely constrain the 
spin-down luminosity of the putative pulsar, the power avail- 
able for particle acceleration is high. In the favoured scenario 
of £ 3 x 10 37 ergs _1 , assuming that the power in injected 
electrons constitutes only a fraction 10~ 2 of the spin-down lu- 
minosity (e.g. Sierpowska-Bartosik & Torres 2008), the power 
in relativistic electron would be larger than our requirement of 
L inj = 2x lO^ergs" 1 . 

Another possibility is that instead of a magnetospheric 
origin, the emission detected by FermifLAT has an Inverse 
Compton origin. In this case the efficiency between the spin- 
down luminosity of the putative pulsar and the observed IC lu- 
minosity is expected to be of 1-10%. Considering the observed 
luminosity of a few times 10 35 erg s _1 , the pulsar would need to 
have a spin-down luminosity of the order of 10 37 erg s _1 . 
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The higher values of spin-down luminosities pose a problem 
for the binary pulsar scenario regarding the shape of the detected 
radio emission. Dhawan et al. (2006) strongly argued that the 
rotating cometary-like tail they detected in VLBA observations 
was related to the shocked pulsar wind contained by the stellar 
wind. Romero et al. (2007) already questioned the ability of the 
stellar wind to contain the pulsar wind from the results of three- 
dimensional hydrodynamical simulations even when consider- 
ing a significantly lower spin-down luminosity of 10 36 ergs _1 . 
The shape of the interaction region between the stellar and pul- 
sar winds is determined by the ratio of wind momentum fluxes 



the flux peak observed at <p — 0.62, the optical depth of yy ab- 
sorption is practically independent of inclination at r 



yy 



Epsr 



(13) 



where M* and v w * are the mass-loss rate and wind velocity 
of the Be star and £psr is the spin-down luminosity of the 
pulsar. For the polar wind, we can take = 10 _8 M o yr _1 
and v w * = 10 3 kms _1 , resulting in a wind momentum flux 
ratio of ?7 W = 0.53(£ , PSR /10 36 erg s" 1 ) (Romero et al. 2007). 
Ignoring orbital motion (which on the other hand would give 
rise to small deviations in the inner region of the interaction, see 
Parkin & Pittard 2008), the opening half-angle of the interaction 
cone would be 



8 = 180° 



1 



(14) 



as seen in analytical models of wind interaction (e.g. 
Antokhin et al. 2004). For a spin-down luminosity of 
10 36 ergs~' the stellar wind would barely dominate the 
pulsar wind (?7 W = 0.53, 8 = 62°). Considering the favoured 
high-altitude location for gamma-ray production (i.e., slot 
gap or outer gap models) or an IC origin of the gamma-ray 
emission, the pulsar spin-down luminosity is of the order of 
Epsr ~ 10 37 ergs _1 . Given this luminosity, the pulsar wind 
would dominate the interaction and would not be contained 
by the stellar wind. The ratio of wind momentum fluxes 
would then be of rj^, 16 and the half-opening angle would 
approach 180°, posing a completely different scenario from 
that usually pictured by wind interaction models (e.g. Dubus 
2006; Sierpowska-Bartosik & Torres 2009; Zdziarski et al. 
2010). A lower value of the beaming correction factor /q in the 
magnetospheric scenario could mean that the pulsar spin-down 
luminosity is as low as £psr - 2 x 10 35 ergs _1 and thus 
contained by the stellar wind. However, given the high energy of 
the cutoff in the Fermi/LAT spectrum of LS I +61 303 we favour 
a high-altitude location for gamma-ray production, leading to 
the higher value of the spin-down luminosity. 



0.2 

for 500 GeV photons. However, the second, wider outburst at 
phases 0.8-1.1, does have an absorption dependency with incli- 
nation, with T yy 1.4 for i = 60° but r yy « 0.3 for i = 15°. 
The diminished absorption at low inclinations leads to a VHE 
flux higher than the observed values, particularly taking into ac- 
count the upper limit measurement at phase 0.8. However, given 
the high uncertainties of the VHE nightly fluxes, the difference 
between the light curves computed for i = 60° and / = 15° is of 
less than l<x. Even though we do not consider it significant given 
the present sensibilities, deeper VHE observations with firm de- 
tections at these phases could help to constrain the inclination of 
the orbit. 



5.4. Magnetic field 

The X-ray/VHE fluxes can be explained using a constant mag- 
netic field of B = 0.22 G. A constraint on the accelerator max- 
imum energy can be put by requiring that the Larmor radius 
of the accelerated electron (r^ = E e /qBc) is contained within 
the accelerator size. This translates to a maximum energy of 
£max = 300 BqR\2 TeV, where Be is the magnetic field in Gauss 
and R12 the accelerator radius in units of 10 12 cm. For the mag- 
netic field we found and considering an accelerator size of the 
order of 10% of the orbital separation, the maximum energy is 
well above the required 10 TeV and radiative losses dominate at 
high energies (see Sec. 3.3). In the binary pulsar modelling ap- 
proach of Dubus (2006), the magnetic flux was considered vari- 
able along the orbit as a consequence of the variation in the dis- 
tance of the wind-shock region to the pulsar. In the shocked wind 
region, the magnetic field values obtained were as high as 10 G. 
Chernyakova et al. (2006) obtain a very similar magnetic field to 
ours, with B = 0.25 G (B = 0.35 G) for the low (high) flux state 
of the source. However, the model is based on an emitting parti- 
cle spectrum very soft at high energies, resulting in a synchrotron 
component peaking in the radio band and an IC component dom- 
inating from soft X-ray to the MeV band. They require another 
hadronic component to explain the VHE emission. The magnetic 
field in this model traces the flux ratio between radio and X-ray 
fluxes instead of ours, which traces the flux ratio between X-ray 
and VHE. On the other hand, in the microquasar approach used 
by Bosch-Ramon et al. (2006b), the magnetic field at the base of 
the jet was taken as 1 G. Most of these approaches use higher 
magnetic fields than the one we found to best reproduce the ob- 
servational features. An explanation for this is likely the need in 
these models to generate both the pedestal and X-ray flux that 
we considered independently. 



5.3. Orbital inclination 

Synchrotron emission will have no dependence on the orbital 
inclination because we are assuming isotropic particle distri- 
bution and a fully disordered magnetic field. For IC emission, 
on the other hand, there is a dependence with the angle 9 be- 
tween the line of sight and the direction of propagation of the 
seed photons (i.e., the connecting line between the star and the 
compact object). Given the orbital configuration of LS 1+61 303 
(see Fig. 1), 9 will not vary significantly with inclination for 
the phases around apastron, at which the source is detected in 
VHE. For yy absorption, on the other hand, the dependency is 
on the angle between the directions of the VHE photon and the 
seed photon integrating along the line of sight out of the system. 
This means that absorption will show weaker orbital variabil- 
ity for low inclinations and stronger for high inclinations. For 



6. Summary and concluding remarks 

We presented a radiation model of a single leptonic population 
that can successfully describe the data obtained through a si- 
multaneous X-ray/VHE campaign on LSI +61 303 performed 
by MAGIC in 2007. The observed X-ray/VHE correlation in- 
dicates both that the emission at these bands may come from a 
single particle population and that energy losses may be domi- 
nated by adiabatic losses. By assuming a constant magnetic field 
and injection along the orbit, we inferred the adiabatic loss rate 
from the X-ray light curve and the time dependent electron maxi- 
mum energy as the balance between acceleration and energy loss 
time scales. We found that a quite efficient accelerator (77 ~ 10) 
is required to obtain the observed VHE spectra because of the 
fast adiabatic cooling. The injected electron spectrum is constant 
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along the orbit and initially taken as a power law with a high- 
energy cutoff, but the Fenm/LAT HE spectrum poses a constraint 
on the hardness of the spectrum at lower energies, which we es- 
timate has to be harder than a e 1.8 for electron energies below 
4 x 10 11 eV. At higher energies, on the other hand, an index of 
a e = 2.1 matches the observed X-ray and VHE photon indices. 
The observed light curves are best reproduced using a magnetic 
field of B — 0.22 G. The general picture shows that the observed 
emission in X-ray and VHE is fully compatible with originating 
in the same parent particle population under dominant adiabatic 
losses. 

The required luminosity budget in injected electrons for the 
computed X-ray and VHE emission is 2 x 10 35 erg/s. Both the 
microquasar and binary pulsar scenarios are able to provide 
these levels of luminosity in accelerated electrons. The GeV 
luminosity of LS I +61 303 can also give us hints on the avail- 
able power: independently of whether the detected emission is 
magnetospheric or IC emission, the injected power should be 
^ 10 37 erg/s. If it is magnetospheric emission, the high-altitude 
favoured scenario of gamma-ray production implies that the 
spin-down luminosity of the pulsar is E — 3 x 10 37 erg/s. This 
high spin-down luminosity poses problems for the binary pul- 
sar scenario because the pulsar wind would not be contained by 
the stellar wind, which would give rise to a completely different 
scenario from the one usually considered. 

The VHE and X-ray light curves of LS5039 and 
LS I +61 303 seem relatively similar, both peaking close to the 
apastron. For LS 5039 the maximum is also close to the inferior 
conjunction of the compact object, whereas this is not the case 
in LSI +61 303. If the light curve maxima are related to apas- 
tron rather than to the observer orientation with respect to the 
system, it may point to a similar physical origin for the X-ray 
and VHE modulation in both sources. On the other hand, the 
adiabatic losses could have a different cause. Recall the very dif- 
ferent stellar winds in both sources. 

We found that the X-ray and the VHE emission are consis- 
tent with a single parent particle population, whereas the HE 
gamma-ray emission detected by Fermi/LAT probably has a dif- 
ferent origin. The orbital anti-correlation and the exponential 
cutoff spectrum suggest that several non-thermal emitting re- 
gions may be present in the system. This could also explain the 
presence of a pedestal emission in X-ray that is compatible with 
constant emission along the orbit. 

We conclude noting that simultaneous observations at the X- 
ray and VHE bands provide a unique diagnostic tool to probe 
the emitter properties in galactic variable sources. Future ob- 
servations, in particular by more sensitive next-generation VHE 
observatories, will allow us to better understand these peculiar 
systems. 
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